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ABSTRACT 

Feedback from outflows driven by active galactic nuclei (AGN) can affect the distribution and 
properties of the gaseous halos of galaxies. We study the hydrodynamics and non-thermal 
emission from the forward outflow shock produced by an AGN-driven outflow. We consider a 
few possible profiles for the halo gas density, self-consistently constrained by the halo mass, 
redshift and the disk baryonic concentration of the galaxy. We show that the outflow velocity 
levels off at ~ 10 3 kms _1 within the scale of the galaxy disk. Typically, the outflow can 
reach the virial radius around the time when the AGN shuts off. We show that the outflows 
are energy-driven, consistently with observations and recent theoretical findings. The outflow 
shock lights up the halos of massive galaxies across a broad wavelength range. For Milky Way 
(MW) mass halos, radio observations by The Jansky Very Large Array (JVLA ) and The Square 
Kilometer Array (SKA) and infrared/optical observations by The James Webb Space Telescope 
(JWST ) and Hubble Space Telescope (HST ) can detect the emission signal of angular size 
~ 8" from galaxies out to redshift z ~ 5. Millimeter observations by The Atacama Large 
Millimeter/submillimeter Array (ALMA) are sensitive to non-thermal emission of angular size 
~ 18" from galaxies at redshift z < 1, while X-ray observations by Chandra, XMM-Newton 
and The Advanced Telescope for High Energy Astrophysics (ATHENA) is limited to local 
galaxies (z < 0.1) with an emission angular size of ~ 2'. Overall, the extended non-thermal 
emission provides a new way of probing the gaseous halos of galaxies at high redshifts. 

Key words: shock waves - galaxies: active - galaxies: haloes - quasars: general - radio 
continuum: general. 


1 INTRODUCTION 

Outflows from ac tive galactic nuc l ei (AGN) regulate b lack 
hole (BH) growth dSilk & ReeJ 199^ DiMatteo et al 20051) and 
may quench star formation JSpringel et al.1120051 : iHopkins et all 
l2008h in galaxies. A great amount of observational evidence 
has demonstrated the presence of AGN-driven out flows, in- 
cluding observations of absorptions in quasars (|Ganguly_etal. 


2011 


2015). 


20071 : [Fu & Stocktonl 120091: iMoe et alj 120091: IVillar-Martm et al 


Arav et alT 20131 : ' Zakamska & Greene! 201 dl T Arav et al. 

_ multiphase outflows in nearby ultraluminous infrared 

galaxies (UL IRGs) (|RuEkg_&Willgux 201 ll: ISturm et al.1 1201 ll : 
ICicone et alj 2014^ Feruglio et alj 2015; Tombesi et al J 20151) and 


quasars dCamian i et alj 1 20151: Gofford et al .1 1 20151) . and post- 
starburst galaxies dTripp et al J|2011|) . The velocity of AGN-driven 
outflows can reach ~ 10 3 kms -1 on galaxy scale, indicating that 
the outflows are likely to propagate into the halos of galaxies while 
the AGN is active. Here we propose to use AGN-driven outflows as 
a probe of the halo gas in galaxies. 

Halo gas has been identified in multiphases (see review 
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by Putman et al. 2 012): cold neutral hydrogen detected as 
high velocity clouds dKalberla et alj 1 20051 : IWestmeier et alj|2005l : 
lOosterloo et al.1 [20071 : Saul et all 20121) . warm gas (T ~ 10 4 5 
K) discovered in deep H a emission line surveys dPutman et al.1 
l2003l : lLehner et al j|2012|), warm-hot gas ( T ~ 10 5-6 K) detected 
in absorption Iwakker & Savage 2009|; Prochaska & Hennawil 


2009i: ISavage et al.ll201 ll: Marasco et alj|20li: F; 


anna et al 


1 JI 2 OI 4 I: 


Ford et all 120141) and hot gas (T ~ 10 K) inferred from X- 


ray observations in emis s ion and absorption Bogdan et aP l2013l : 
iMiller & Bregmanl 120131 : iBogdan et all 12015 ). The presence of 
warm-hot and hot halo gas, extending out to the virial radius, is 
of particular interest since the hot gas is postulated to host a sig - 
nificant fraction of baryons in the galaxy dKaufmann et al Jl2006l ). 
However, the detailed properties and the origin of the extended and 
diffuse hot halo gas remain uncertain si nce th ere is_ littl e evi dence 
for its existence around spiral galaxies dPutman et al J|20 1 2 |) . The 
detection of halo gas out to virial radius scale is difficult and the 
extent to which the outflows impact the properties of the halo gas 
remains uncertain. Therefore, it is important to study the interac¬ 
tion between AGN-driven outflows and surrounding gas on differ¬ 
ent scales as a probe of the properties of the diffuse hot halo gas 
and the effectiveness of the feedback mechanism. 
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2 X. Wang and A. Loeb 


In galaxies with a weaker AGN where the energetics of 
AGN activity and star formation are comparable, it remains un- 
cl ear w hether outflows are dominated by AGN or supernovae (SN) 
dHonkins et aljEoT^ l. In this paper, we focus on AGN-driven out¬ 
flows. First, our model assumes spherical symmetry, which is more 
justified for AGN-driven outflows since they are launched at the 
center of the galaxy whereas SN-driven outflows are distributed 
throughout the entire disk. More importantly, as shown later in 
the paper, the strongest emission signal comes from more massive 
galaxies where AGN feedback is thought to dominate. 

Previous work on the dynamics of and non-thermal emission 
from galactic outflows has made simple assumptions about the total 
gravitational mas s and the gaseous environment in wh i ch the out¬ 
flow p ropagates l Furlanetto_&Loeb 200l|; King 20031; King et all 


l201ll : iFaucher-Giguere & Ouataertf 20121 ; iNims et aTfl2015h . and 
limited the evolution of the outflows to galactic disk scales 
1 Jiangegak 2010|^ Faucher-Giguere & Ouataerd 120121 ; INims et aTI 
120151: iHopkins et al.ll2015ir In this paper, we explore different gas 
density profiles in galaxy halos and examine the non-thermal emis¬ 
sion from the forward shock plowing into the ambient medium in 
details. We predict the multiwavelength spectrum and detectabil¬ 
ity of the non-thermal emission and discuss how the outflow shock 
and halo gas affect each other. We propose a new way to probe 
the gaseous halo using the non-thermal emission from the outflow 
shocks as they travel through the ambient medium in the galaxy 
and halo. 

The paper is organized as follows. In § [2] we describe our 
model for the halo and gas distribution. In § [3] we analyze the hy¬ 
drodynamics of AGN-driven outflows. In §[4] we calculate the non- 
thermal emissions from shocks produced by outflows. In § [5] we 
show numerical results for representative cases and discuss physi¬ 
cal significance. Finally in § [6] we summarize our results and dis¬ 
cuss their implications. 


2 MODEL DESCRIPTION 

We approximate the galaxy and halo as spherically symmetric. The 
environment into which the outflow propagates is decribed below. 
Here we di scuss properties of spherical outflows dr iven by fast nu¬ 
clear wind jjiang et all20ld ; lKing & Poundsll2QI5l) . The predicted 
radio emission from outflow shocks as discussed in § [4] is fainter 
than the radio synchrotron e mission from relativistic jets in a small 
subset of all active galaxies l iHeckman & Bestll20l4 ). 


2.1 Mass profile of host galaxy 


We assume that the density distribution of the galaxy in which 
the o utfl ow is ini tially embedded follows the NFW profile 
dNavarro et al 11993) : 


Pnfw (-^) Po(l z) 




flm(z) C N x(l + C N X ) 2 


(1) 


where po = 3 Hq/8tvG is the critical density today, Ho is 
the Hubble constant today, G is the gravitational constant, x = 
R/R v ii, c N is the concentration parameter which is roughly given 
by: c N w 25(1 + a) -1 , Sl m = 0.3. (5 C is given by S c = 
A c ci,/[3(ln(l + c N ) - c N /(l + c N ))}, where A c « 187r 2 . 
flm(t) can be expressed as fi m (z) = fi m ( 1 + 2 ) 3 /[fi m ( 1 + 
z) 3 + Ha], where Ha = 0.7. R v i r is the virial radius, writ¬ 
ten as R vit = 0.78 h~ 2/3 [fi m A c /187r 2 fi m (2)]~ 1/3 M 3/3 /{l + 
zj 10) kpc, where h = (Tfo/lOOkms -1 ) is the Hubble parameter 


and A/haio = 10 8 Mh,8Mg is the halo mass. We obtain the total 
mass of the galaxy and dark matter halo within a radius of R by 
f 4-7T R 2 p NFW (R) dR , which gives: 


Mqm(R) =Po(1 + z) 


3 H m ScRvh 


Hm (z) 


ln(l + c N x) - 


N 

c„x 


1 + c N x 


( 2 ) 


We es timate the BH mass M, self-consistently by the following 
steps dGuillochon & Loebll2015h . First we o btain the total stell ar 
mass in the galaxy M* determined by AThaio JMoster et alj201(il) : 


Mi, = 0 - 


(AThaio/ATi 


[l + (M halo /M 1 ) /3 ] 


(71 — n)/P ’ 


(3) 


where log(M* ! o/MQ) = 10.864, log(Mi/Af 0 ) = 10.456, 71 = 
7.17, 72 = 0.201 and /3 = 0.557. There is no specific bulge mass 
A/buige for a given halo mass dKormendv & I hill 201 ill . Numerical 
simulation dBluck et alj|2014h suggests that the bulge-to-total stel¬ 
lar mass ratio B/T =Mbui ge /Af* is roughly uniformly dist ributed 
from 0 to 1. This ratio for the MW is ~ 0.15 dLicquia & Newman] 
l2014lt . Additionally, Fisher & Drory (2011) suggest that ~ 25% of 
all local stellar mass is in bulges and elliptical galaxies. We then 
adopt a particular value of B/T ratio to be ~ 30% in our calcula¬ 
tion and multiply this value by AT* to get ATbuige to illustrate some 
examples. There is likely to be only ellipticals in high mass halos, 
so it is justified to take a fixed B/T ratio for these systems. We also 
verify that modifying B/T ratio only results in a difference within a 
factor of 4. This variation can be cancelled out by the uncertainty 
in the fraction of AGN’s luminosity injected into the medium as 
discussed later in the pap er. Finally, we obtain the BH mass M. by 
dMcConnell & Mall2013l> : 


log (M./Mq) = 8.46 + 1.05 log 


ATbulge 

lO n Af 0 


(4) 


The underestimation of the M, correlations co uld be an issue for 
the most massive BHs dKormendv & H0II2OI3I) but should not af¬ 
fect our results on the emission from the outflow shocks in the more 
common galaxies with M. 10 9 AT 0 . 


2.2 Gas density profile 


We assume that the gas takes up a fraction / g of the total mass 
of the dark matter in a galaxy. We adopt_a cosmic mean baryon 
fraction, which is f s ~ 16% dHinshaw et aLlfioiBl) . A fraction of 
the baryons /a is concentrated in the disk of the galaxy, and the 
d isk r adiu s Rusk is taken to be ~ 4% of the virial radius R v i r 
dShibuva et alj[2015h . 

Our first prescription for the gas density distribution is a bro¬ 
ken power-law profile, given by: 


ppi(R) 


C d R~ a (R ^ .Rdisk) 

Ch R^ 13 (77disk < R ^ Rv ir) 


(5) 


where a and f3 are the power-law indices in the disk and halo com¬ 
ponent, respectively. We assume an isothermal sphere for the gas 
within the disk component and fix a = 2.0 in our calculation. The 
constants in the density profile C d and Ch can be constrained by 
the baryon mass budget in the disk component and in total. Con¬ 
sequently, p is soley dependent on f d . The constraint on /3 by f d 
is shown in Fig[T] where we find that when f d ~ 0.25, j3 ~ 3, 
indicating that the gas in the halo approximately follows the NFW 
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Probing gaseous galaxy halos 3 



Figure 1. Power-law index f) of the halo gas density profile as a function 
of the baryon fraction of the halo (1 — /<j). The dashed lines correspond 
to values of f,\ = 0.1, 0.25 and 0.5, which we have taken into numerical 
calculation in the following sections. 


profile. From the broken power-law density profile, we estimate the 
gas number density at 50—100 kpc to be 1 0~ 5 —10~ 4 cm~ 3 , which 
is consistent with numerical simulatio ns iSokoIowska et al.ll2015tl 
and observations iBogdan et al J|2Q 1 5l) of the hot halo gas distribu¬ 
tion. 

The second profile we consider for the halo gas density 
distr ibution is analogous to that of galaxy clusters, written as 
iPatei & Loehll2015l) : 


Pclu(f?) 


r f s A 


(R/s) 2T ~ 2 

R/r s [1 + (s/r s )(f?/s) r ] 2 ’ 


(6) 


where A = po5 c is the scale parameter, s = R v i r , r s = s/cn is the 
scale radius and T is the jump ratio. The density profile recovers to 
a scaled NFW profile for T = 1. 


3 HYDRODYNAMICS 

We assume spherical symmetry for the outflow and the ambi¬ 
ent medium. Fast wind with velocity ~ 0.1c is injected into the 
medium, as inf erred from observations of broad absorption lines in 
quasars dArav et al ■ 2011. The wind drives an outer forward shock 
into the ambiet medium accelerating the swept-up material and an 
inner reverse shock into the wind decelerating itself, separated by a 
contact discontinuity iKing & Poundslfioil . 

The equation of motion of the shell is gi ven _ by 

dFurlanetto & Loebl 1200 ll : IFaucher-Giguere & Ouataertl l2012h : 

d 2 R s _ 4ivR 2 GMtot v s dMs 

^-^ir {Ft ~ Fo) '~Rr~w B ^r' (7) 

where G is the gravitational constant, and R s , v s and M s are the 
radius, velocity and mass of the swept-up shell, respectively. Mtot 
is the total gravitational mass inside R s that impedes the expansion 
of the wind bubble, written as Mtot = Mom + M ga i + M. + 
M s / 2, composed of the mass of dark matter Mdm, galaxy M ga j, 
the central BH M,, and the self-gravity of the shell. The shell mass, 
M b , satisfies, 

dMs . 

—^-=4np g R s Vs, (8) 


where p g is the ambient gas density profile in the galaxy. 

Hydrostatic equilibrium gives the temperature in the ambient 
medium To: 

dTo _ GMtotmp T 0 dn g 

Hr ~ kR? n//~dR ’ ( ’ 

where m p is the proton mass, k is the Boltzman constant and n g is 
the number density profile of the ambient gas. At virial radius R v i r . 
To reaches virial temperature T v i r = pm p v 2 /2k where p = 0.5 
is the mean molecular weight of fully ionized gas and v c is the 
circular velocity, given by v c = (GMh a io/T v ir) 1//2 . The ambient 
thermal pressure is given by Pq = n g kTo. 

The thermal pressure in the shocked wind Pt declines due to 
radiative energy losses and work done on the ambient gas by the 
expansion, at a rate: 


dPt 

dt 


A 

2ivR3 



( 10 ) 


where A is the heating and cooling function, composed of energy 
injection from the central source and different physical cooling pro¬ 
cesses in the shocked wind region: 


A — Tin Tff Tic T S y n T p . (11) 


Energy is continuously injected into the shocked wind during the 
quasar’s lifetime, taken to be the e-folding time TEdd ~ 4.5 x 
10 7 yrs d Martini & Weinbergll200ll) . with a rate of Ti n , which is 
assumed to be a fraction of the AGN’s bolometric lumi nosity 
.finTAGN- Observations infer /i n to be ~ 1% — 5% dArav et al .1 
120131 : ICicone et al .Il2014h and we adopt /i n = 5% in our calcula¬ 
tion. We assume that Tagn is a fraction / AGN of the Eddington 
luminosity TEdd = 1-38 x 1 0 38 (M. /Mp,) e rgs^ 1 , and adopt 
/agn = 0.5 in our calculation dShen et alj2009l) . 

The last four terms in the right hand side of Eqn. QT] ac¬ 
count for radiative cooling. Tff is the radiative cooling rate via 
free-free emission in the shocked wind. Tic decribes cooling via 
inverse Compton (IC) scattering off photons in the quasar’s radi¬ 
ation field and the cosmic microwave background (CMB). T syn 
represents synchrotron cooling rate. T p refers to the cooling of 
protons through Coulomb collisions with the electrons. The cool¬ 
ing rate can be expressed as pEt/t c , where Et = 2ivR^Pt is 
the thermal energy in the shocked wind and t. c is the timescale 
corresponding to different cooling p rocesses. The total emissiv - 
ity of free-free emission is given by dRvbicki & Liehtmanlll979t) : 
£fi = 1.4 x 10~ 22 T^o n 2 0 g B , where g B is the Gaunt fac¬ 
tor, T e ,io = (T e /10 10 K) and n e ,o = (n e /lcrrC 3 ) are the 
electron temperature and number density, respectively. The cor¬ 
responding cooling timescale is ts = |fcT e /eff = 4.69 x 
10 8 T^ 2 n~g g B 1 yr. The IC cooling time of relativistic elec- 
tro ns of energ y E e in soft photon radiation field can be written 
as ( iKing & Pounds! lit) 1 5i) : tic = 3m 2 c 3 /8iraTU p hE e , where 
or is the Thomson scattering scross section and f/ p h is the en¬ 
ergy density of soft photons, including AGN photons with en¬ 
ergy density Tagn = Tagn/47tJ? 2 c and CMB photons with 
energy density Tcmb ~ 4.2 x 10 _13 (1 + z) 4 erg cm -3 . Here 
we consider the most efficient IC cooling limit and thus leave 
out non-relativistic elect rons, o f which the IC coolin g time can 
be significantly longer ( IFaucher-Giguere & Ouataertl 201211 . We 
obtain the temperature in the shocked wind by the Rankine- 
Hugoniot jump condition T e ~ 3pm, p v 2 n /16k. The synchrotron 
cooling timescale is given by f syn = 1.6 x 10 12 BZg T~/ 0 yr, 
where B -6 = (T?/10~ 6 G). If two- tempe rature plas ma ef - 

fect is taken into account iFaucher-Giguere & Ouataertl 1201 2ll . 
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4 X. Wang and A. Loeb 


then the proton cooling timescale f p can be expressed as: 
t p ft: 1.4 x 10 9 i^ ikpc LXgn ,46 v s,3 5 y 1 - where «s,3 = 

(tts/10 3 kms -1 ) and Lagn,45 = (Lagn /10 43 ergs -1 ). 


4 NON-THERMAL EMISSION 

Next we discuss the non-thermal emissionjrom th e out flow shock 
as it propagates in the ambient medium iNims et al.l2015t) . 


4.1 Synchrotron emission 

As the forward shock plows through the ambient medium super¬ 
sonically, a broken power-law distribution of non-thermal elec¬ 
trons N(yy) drf oc -y~ p (1 + 7/7b)~ 1 is generated via Fermi ac¬ 
celeration in the shock to produce non-thermal emission, where 
p is the power-law index. 7b is the break Lorentz factor, which 
is obtained by equating the dynamical timescale ~ Rs/v s and 
the cooling timescale 3m e c/4(C/B + Uagn + Fcmb )ot7- This 
gives 7b = 3m e cv s /AaTR a (UB + Uagn + Ucmb), where m e 
is the electron mass, <tt is the Thomson scattering cross section 
and Ub = B 2 /8tt is the energy density of the magnetic field. We 
assume that the total non-thermal luminosity is a fraction of the ki¬ 
netic energy of the swept-up material, written as L n t = e n tLkin ~ 
ie n t M s v 2 . We calibrate the magnetic field energy density as a frac¬ 
tion £ b of the thermal energy behind the shock in what follows 
supernova (SN) remnants l lChevalieill 19981) . giving: 

U B = £, B nkT . (12) 


Observations of radio emitting bubbles from a radio-quiet quasar 
imply p ~ 2 OHarrison et al .1120151 ). By fitting the radio flux from 
bubbles at ~ 10 kpc, we obtain e n t ~ 5%. Coefficients £ B can 
be estimated from observations of late-time radio e missio n f rom 
relati vistic jets associated with tidal disruption events dBower et al J 
|201 3l) . synchtrotr on emission from sho cks between jet and circum- 
nuclear medium l iMetzger et al.l | 2012h as well as from an anal¬ 


ogy with SN remnants iChevalien 19981) . These observations imply 


Cb ~ 0.1. 

Finally, we calculat e the synchrotron emission following the 
stand ard formula from dPacholczvkl Il97d : iRvbicki & Lightmanl 
Il979h . The emission and absorption coefficients are given by: 


•syn 

Ju 



F(x)N(-y)d'y , 


(13) 


= —C2 B- 


/ , ' ymax 2 d 

[JV(7)1 

A min ' 

. 7 2 . 


F{x) d.'y , (14) 


where ci = v / 2e 3 /47rm e c 2 , C 2 = \/2e 3 / 8-ivm 2 c 2 , F(x) = 
x f°° -ff 5 / 3 (£) and K 5 / 3 (x) is the modified Bessel func¬ 
tion of 5/3 order. The maximum Lorentz factor 7 ma x is given 
by the tighte r constraint of equaling th e acceleration timescale 
^b-ccRbc/v 2 dBlandford & Eichlerl Il987h to either dynamical or 
cooling timescale, where ( a cc ~ 1 and Rb = 7 m e c 2 /eB is the 
Larmor radius. We plot 7 max in unit of 10' as a function of outflow 


shock radius R B for Mh a io = 10 12 Mg, /a = 0.25 and 2 = 1.0 
as a representative example, shown in Fig|2] 7 ma x varies within a 
factor of ~ 5 as a result of simultaneously decreasing v s and soft 
photon energy density with increasing R s . We take the minimum 
Lorentz factor 7 m in ~ 1 in our calculation. The synchrotron emis¬ 
sion peaks at a frequency of u syn = 4.2 x 10 14 73-6 7 ? Hz, where 
77 = (7/10 7 ). 



Figure 2. The maximum Lorentz factor of non-thermal electrons 7max 
in unit of 10 7 as a function of outflow shock radius. We fix Mb a i 0 = 
10 12 Mq, j\t = 0.25 and z = 1.0 as a representative example. 


4.2 Inverse Compton scattering 


The soft photons includes those from the accretion disk and CMB. 
The energy density of the AGN radiation field is (7 agn ~ 2.8 x 
10 -1 °Lagn, 45 R s ^ pc erg cm 3 . The CMB photons have an en¬ 
ergy density of U CMB oc (1 + 2) 4 , which manifests them selves as a 
dominant source of IC scattering at high-redshift dCelotti & Fabian! 
1 20041) . The spectral ener gy distribution of quasars can be con¬ 
strained by observati ons J E1 v i s et ali i 19941: iMarconi et all |2004 
IScott & Stewartl2014l). For_simplicity, we approximate it as a black 
body spectrum l lto et al.| 120151) . We model the CMB photons as 
a black body with a spectrum peak frequency of tz CMB fs 1.6 x 
10 n (l + 2) Hz. The peak of IC scattering of CMB photons takes 
place at a frequency of u lc « 7 2 ^ CMB = 1.6 x 10 25 77(1 + 2) Hz. 
The differential ra te to produce high-energy photons wit h energy 
em e c 2 is given by jjonesll 1 9+ilConni & Blandfordlll990l) : 

Q(e) = J de 0 n{t 0) J d'yNt'y) K(e, 7, e 0 ) , (15) 

where eo m e c 2 is the soft photon energy, 7 m e c 2 is the electron en¬ 
ergy and n(e 0 ) is the number density of soft photons. K(e, 7 + 0 ) 
is the Compton kernel, expressed as: 


p7T7*^ n 

K(e, 7 , eo) = [2k In k + (1 + 2k) (1 - k ) 

7 £ o 


+ 


(4ep7K) 2 

2(1 + 4eo7«0 


(16) 


(1 - «)] , 


where k = e/[4eo7(7 — e)]. The emission coefficient of IC scat¬ 
tering can be obtained by: 


•IC 


= ^ 


(17) 


where h is the Planck constant. 


5 NUMERICAL RESULTS 

In Figures[3}j6] we show the dependence of outflow hydrodynamics 
solutions and emissions on /a, Mhaio, z and density profile formu¬ 
lation. Since the gas distribution in the intergalactic medium (IGM) 
is uncertain, we restrict our calculation to halo scale within R v n. 


© 2015 RAS, MNRAS 000 .mfl 2 l 













































Probing gaseous galaxy halos 5 


As shown in panel a in Figures [3}j6] we find that the swept- 
up shell decelerates quickly to a roughly constant velocity of ~ 
10 3 kms -1 in the disk. As it propagates outside the galaxy into 
the halo, the shell accelerates somewhat as a result of the tenu¬ 
ously distributed halo gas. The evolution of the shell velocity is 
consistent with a self-similar solution, where the shell radius is as¬ 
sumed to follow R s oc t s and v s oc t lS_1 . We express the gas 
power-law density profile generally as p oc R 1 . For 7 < 3, we 
obtain M a oc R 3 T In the energy-conserving limit, we assume 
that ~ 50% of the injected energy goes to the kinetic energy of the 
swept-up material, Lint = M s v a , and so we have 5 = 3/(5 — 7 ). 
For power-law index a = 2 in our model, 5=1 and thus v s ap¬ 
proaches a constant in the disk. We can also verify that for halo 
component power-law index /3, the outflow accelerates as /3 > 2. 
The acceleration stops as the quasar shuts off and the thermal en¬ 
ergy in the shocked wind E t drives the expansion of the shell after¬ 
wards. At this point, the outflow reaches the edge of the dark matter 
halo and is likely to continue to propagate into the IGM. 

Panels b and c in Figures [3j|6] show the radio flux as a func¬ 
tion of shock radius and time, respectively. We scale the time 
to the Flubble time tn, which is given by tn = 1/H(z) = 
Hq 1 [fi m (l + z) 3 + 12a] 1 ^ 2 . The chance of finding a galaxy 
with a given flux is t/tn. We find that for z ~ 1, about a few per¬ 
cent of the galaxy halos embed outflows reaching R v i r . We also cal¬ 
ibrate the angular diameter of the outflow shock, given by 7t! a / Da, 
where Da is the angular distance. 

We show snapshots of non-thermal emission taken at two 
milestones in panels e and/. At the edge of the galaxy disk, the en¬ 
ergy injection from the central source has an age of ~ 10 ' yrs. At 
the virial radius, snapshots are taken at the dead quasar remnants 
with outflow approaching the edge of the dark matter halo on a 
timescale of ~ 10 s yrs, which indicates that this population should 
be ~ 10 times more abundant. At this point, the outflow no longer 
overlaps with the galaxy and there is no galaxy-bubble interactions. 
We find that the outflows can reach the edge of the halo around the 
end of quasar’s lifetime. This feature indicates that AGN-driven 
outflows are most abundant during their passage through their host 
galaxy halo. 

We summarize the detectability of this extended non-thermal 
emission in Table [Tj 


5.1 Dependence on parameters 

5. 1.1 Disk mass fraction 

For a halo of mass ALhaio = 10 12 Mq at z = 1.0, we choose 
three re presentative values of /a as motivated by observations 
dCourtois et al.ll2ol5l) . We find that the shell velocity is not sen¬ 
sitive to /a. The outflow reaches the edge of the halo around 
the time the energy injection discontinues. With a velocity of 
~ 500 — 10 3 kms -1 , the outflow is likely to propagate into the 
IGM. The non-thermal radio flux at 1 GHz remains at ~ 0.1 mjy 
within the disk, independent of /a. As the shell propagates into the 
halo, the non-thermal emission diminishes quicker in halos with 
higher /a as a result of more tenuous halo gas. For /a = 0.5, 
the radio emission is ~ 100 times fainter than the other two cases 
and drops below the detection limit of JVLA and SKA before the 
outflow reaches R v n. Observationally, we can distinguish galaxies 
with high disk baryonic concentrations by the faint emission from 
their outflows propagating in the halos. 


5.1.2 Halo mass 

We examine ALhaio of 10 11 Mq, 10 12 Mq and 10 13 Mq, covering 
the full range from mid to high mass halos. In lower mass halos, 
the energy input into embedded outflows is much lower due to the 
self-consistent scaling relation between M . and Mhaio- The out¬ 
flow shock decelerates quicker and may not propagate farther out¬ 
side the galactic disks. The short lifetime of outflows in low mass 
galaxy halos makes them less abundant. Therefore, it would be ob¬ 
servationally challenging to identify outflows from low mass halos 
in terms of both emission intensity and recurrence rate. At 3 ~ 1, 
the emission is only detectable in radio band on galaxy scale with a 
flux ~ 10 p5y. High mass galaxies produce AGN photons of higher 
energy density, making the detection more promising. 


5.1.3 Redshift 

The hydrodynamics of outflows is insensitive to 2. Consequently, 
outflows reach the edge of its host galaxy and halo at similar veloci¬ 
ties for different redshifts. At low redshift z ~ 0.1, the non-thermal 
emission is detectable in multiwavelength from radio to X-ray. For 
high-redshift galaxies at z = 5, the non-thermal emission is dom¬ 
inated by IC scattering off CMB photons. The emission remains 
observable in the radio, infrared and optical bands on halo scale. 


5.1.4 Gas density profile 

We compare the broken power-law profile to the gas density profile 
of galaxy clusters. We find that the outflow velocity and emission 
indistinguishable for these gas density profiles. However, outflows 
can not reach the edge of the halo for galaxy clusters, excluding 
them from halo scale observations in these systems. 


5.2 Energy or momentum conserving outflow 


Another important dynamics issue is whether the outflow is mo¬ 
mentum or energy conserving. In the momentum-driven regime, 
thermal energy in the shocked wind region is efficiently radiated 
away, while in energy-driven outflows, such radiative losses are in¬ 
significant. We compare the timescale of the most efficient radiative 
cooling processes discussed in §3 in the shocked wind, t coo i, with 
the dynamical timescale of the outflow, given by fdyn = R s /v s , as 
shown in Fig. [7] 

The plot shows t coo \/ta yu for several representative cases and 
indicates that for some cases the outflow starts propagating as par¬ 
tially momentum-driven. Once the shell reaches ~ 100 pc, the 
partially momentum-driven regime breaks down and the shocked 
wind region no longer cools rapidly. At larger radii, the soft pho¬ 
ton energy density is dominated by CMB photons and f CO oi/fdyn 
decreases consequently. However, the energy conserving nature re¬ 
mains unchanged at larger radii, w hich is in agreement with re¬ 
cent^observationsjTbmbesGlt alj2015h and theoretical calculations 


jFaucher-Giguere & Ouataertll20121 : Zubovas & Kindl2()T3) . 

These results suggest that most of the wind kinetic energy is 
converted to the kinetic energy of the outflow, giving P 2 /M s ~ 
P 2 a d/AL n , where P ra d = Lagn/c is the momentum flux of 
AGN’s radiation field an d Mj n is the mass inject ion rate of the wind 
from the central source iZubovas & Kingll2012h . We can write the 
momentum flux of the outflow normalized to AGN’s radiation as 
P/Prnd ~ Vi n /v a . This relation is illustrated in panel d of Figures 
3-6. 
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log(i? B /pc) log (t/t u ) 


log(Rjpc) 



log(^/Hz) 


Figure 3. Dependence of outflow hydrodynamics and emission on baryon fraction in the disk /d- We fix Mh a io — 10 12 Mq and z = 1.0. Panel a and b 
show the shell velocity and radio synchrotron flux at 1 GHz as a function of radius. The dotted and dashed vertical lines mark the position of -Rdisk an< l ^vir> 
respectively. The upper x-axis of panel b marks the angular diameter of the outflow shock. Panel c shows the radio synchrotron flux as a function of time. The 
dashed vertical line corresponds to the point when the AGN shuts off. Time is scaled to the Hubble time £h on the upper x-axis. Panel d demonstrates the 
momentum flux boost of the shell. The solid lines represent the numerical result while the dashed lines correspond to predictions in the energy-driven regime. 
Panel e and / illustrate snapshots of non-thermal emission power and flux at .Rdisk and R v i r , respectively. The solid, dashed and dotted lines correspond to 
synchrotron emission, IC scattering of accretion disk photons and CMB photons, respectively. 
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Figure 4. Dependence of outflow hydrodynamics and emission on halo mass Mh a io- We fix /d = 0.25 and 2 : = 1.0. The configuration and physical 
significance of the subplots are the same as Fig. [3] The dotted vertical lines marks the position of -Rdisk f° r the three cases in panel a and b. 
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Figure 5. Dependence of outflow hydrodynamics and emission on redshift z. We fix Afhaio — 1O 12 M0 and /d = 0.25. The configuration and physical 
significance of the subplots are the same as Fig. [3] The dotted vertical lines marks the position of i^disk f° r the three cases in panel a and b. The right-hand 
y-axis of panel e and / is scaled to a distance of 10 Gpc. 
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Figure 6. Dependence of outflow hydrodynamics and emission on gas density profile of galaxy clusters. We fix Mh a io — 10 12 Mq and 2 = 1.0. We compare 
galaxy cluster gas density profile with the broken power-law profile (/^ = 0.25). The configuration and physical significance of the subplots are the same as 
Fig.DO 
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II 

O 

h- 1 


2 = 1.0 


2 = 5.0 


Telescopes 

Fv (T?disk) » Fi/(Rv ir) 

(mJy) 

detectability 

Fv (7?disk) » F u (Rv ir) 

(mJy) 

detectability 

Fv{Rdisk) » Fv(R-v ir) 
(mJy) 

detectability 

JVLA 

300 ; 0.8 

Yes ; Yes 

1.0; 5 x 10“ 3 

Yes ; Yes 

4 x 10“ 3 ; 6 x 10 -4 

Yes ; Marginal 

SKA 

300 ; 0.8 

Yes ; Yes 

1.0; 5 x 10 -3 

Yes ; Yes 

4 x 10“ 3 ; 6 x 10~ 4 

Yes ; Marginal 

ALMA 

0.5 ;4 X 10“ 3 

Yes ; Marginal 

5 x 10" 3 ; 7 x 10“ 5 

Marginal ; No 

4 x 10“ 5 ; 4 x 10~ 4 

No; No 

JWST 

3 x 10“ 4 ; 5 x 10“ 6 * 

Yes ; Marginal 

5 x10“ 5 ;2 x10 -6 * 

Yes ; No 

4 x 10“ 8 ; 2 x 10~ 5 

No ; Yes 

HST 

2 x 20“ 4 ; 5 x 10“ 6 * 

Yes ; No 

3 x10“ 5 ;2 x10“ 6 * 

Marginal ; No 

4 x 10“ 8 ; 2 x 10~ 5 

No ; Marginal 


t^Fj/ (7?disk) 5 oFi/^R-y ir) 
(erg cm -2 s — - 1 ) 

detectability 

vFv (Tt^disk) > ^(tfvir) 
(erg cm -2 s — - 1 ) 

detectability 

vFv(R disk ) , oFvi^R-v ir) 

(erg cm -2 s — - 1 ) 

detecatbility 

XMM-Newton 

10 _io ; 2 X 10 _io 

Marginal ; Marginal 

5 x 10 _la ; 4 x 10“ 1Y 

No ;No 

7 x 10 _2U ; 10~ 1Y 

No; No 

ATHENA 

10“ 16 ; 2 X 10“ 16 

Yes ; Yes 

5 x 10“ 19 ; 4 x 10- 17 

No ; Marginal 

7 x 10“ 2 ° ; 10- 17 

No; No 

Chandra 

2 x 10- 17 ; 7 x 10“ 16 

No ; Marginal 

5 x 10“ 19 ; 7 x 10- 17 

No ;No 

6 x 10 “ 2 ° ; 10- 17 

No; No 

NuSTAR 

2 x 10- 17 ; 7 x 10 -16 

No ;No 

5 x 10 -19 ; 7 x 10 -17 

No; No 

6 x 10 _2 ° ; lO” 17 

No; No 


Table 1. Detectability of non-thermal emission from AGN-driven outflow shock. 

Note: We choose a i Q = 10 12 Mq and = 0.25 as a representative example for a galaxy halo. For radio, mm/sub-mm, infrared and optical 
observations, we provide values of Fv(Rdisk) and F U (R V i r ), which correspond to non-thermal flux at the edge of the disk and halo respectively, in 
unit of mJy. For X-ray observation, we present isF u (.Rdisk) and vFv(Rv i r ), in unit of erg cm -2 s —1 . 

The telescope detection limits are as follows: 

- The Jansky Very Large Array ( JVLA ): ~ 1/zJy for 1 a detection and 12h integration time at most bands JNRAQ|20141) . 

- The Square Kilometer Array (SKA-MID): ~ 0.7/iJy RMS sensitivity for lOh integration time ( IPrandoni & Sevmouil20l4) . 

- The Atacama Large Millimeter/submillimeter Array (ALMA): At observating frequency 345 GHz, the sensitivity ~ 8.7 ply for lOh integration 
time is calculated by the ALMA Sensitivity Calculator (ASC) (https://almascience.eso.org/proposing/sensitivity-calculator). 

- The James Webb Spac e Telescope (J WST): sensitivity ~ 10 nJy for wavelength 1 — 3 pm and ~ 30nJy for wavelength 4 — 5 pm for 10cr detection 
and 10 4 s integration time JSTScll2013l) . 

- Hubble Space Telescope (HST): sensitivity ~ 40 — 50 nJy for wavelength 0.6 — 1.5 pm for 10cr detection and 10 4 s integration time iSTScJ 

12013 I) . 

- Chandra: sensitivity of high resolu tion camera (HRC) ~ 9 X 10 16 erg cm" 2 s" 1 covering energy range 0.08 — 10 keV for 3cr detection and 
3 X 10 5 s integration time ICXdllO i'll) . 

- XMM-Newton: ~ 3.1 X 10~ 16 erg cm -2 s _1 in 0.5 — 2.0 keV band jHasinger et all200lh . 

- Advanced Teles cope for High Energy Astrophysics (ATHENA): ~ 4 x 10 —17 ergcm -2 s —1 in 0.5 — 2 keV band in a 10 6 s deep field 
iBarcons et alkoid . 

- Nuclear Spectro scopic Telescope Array (NuStar): ~ 2 X 10~” 15 erg ran 2 s' 1 in 0 — 10 keV band for 3<x detection and 10® s integration time 
l H arrison et al. l2013il . 

* The emission may be contaminated by scattered quasar light (see §[6}. 


6 CONCLUSIONS AND DISCUSSION 


We study the hydrodynamics of AGN-driven outflows out to galac¬ 
tic halo scales and the resulting non-thermal emission from the fast 
forward outflow shock propagating into the ambient medium. We 
have found that the outflow decelerates rapidly to a nearly constant 
velocity of ~ 10 3 kms -1 within the galaxy disk and accelerates 
once it enters the halo until the central BH shuts off. Around this 
time, the outflow can reach the edge of the halo. We have verified 
that the outflow is energy- conserving on l arge r adii, consistently 
with recent observations (Tombesi^ydJ 20j_5|) and theoretical 
predi ctions l lFauc her-Gigu tire~& Ouataertl 120121 : IZubovas & Kind 
US). The predicted non-thermal emission from outflow shocks 
in MW mass halos up to a redshift z of 5 is detectable over a 
broad range of wavelengths. At z ~ 0.1, the 2' angular scale 
emission is detectable by JVLA and SKA in radio band, ALMA in 
mm/sub-mm band, JWST and HST in optical and infrared bands, 
marginally detectable in X-ray band by Chandra, XMM-Newton 
and ATHENA. At 2 ~ 1, the signal remains observable in radio 
band and marginally detectable in infrared and optical bands with 
an angular scale of ~ 18”. The detection is promising even at high 
redshifts (z ~ 5) in the radio, infrared and optical bands with an 


angular scale of ~ 8”. For lower mass halos the detection should 
limit within the local Universe. 

We find that the detailed gas distributions do not significantly 
affect the hydrodynamics of the outflow while the halo mass plays 
a more important role in regulating the outflow dynamics. We 
show a near universality of the non-thermal emission within the 
galaxy disk for different gas distributions of galaxies with same 
halo masses, which breaks down on the halo scale as a result of 
distinct density profile for tenuous halo gas. The halo mass deter¬ 
mines the intensity of the emission since the BH mass is scaled 
self-consistently with halo mass. Consequently, non-thermal emis¬ 
sion from outflows embedded in low mass halos is ~ 1 — 3 orders 
of magnitude fainter than that in MW mass halos. We conclude that 
the halo mass is the dominant factor in regulating the dynamics and 
emission of the outflow. In order to distinguish between different 
gas density distributions, halo scale observations are required. 

The predicted non-thermal emission should be an observa¬ 
tional signature of the existence of extended gas in galaxy halos in 
a wide range of redshifts. With this method, one can probe the evo¬ 
lution of gaseous halos at early cosmic times. Thermal X-ray emis¬ 
sion from free-free cooling at the forward wind shock was proposed 
to be an observational signature of kpc-scale outflows powered by 
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Figure 7. Ratio of radiative cooling timescale in the shocked wind region to 
outflow’s dynamical time. The default values of the parameters are: J\] = 
0.25, Mh a i 0 = 10 12 Mq and z = 1.0. Each line represents a specific 
parameter modified from its default value while the other parameters are 
fixed at the default values. The dotted line separates the momentum and 
energy conserving regimes. 


AGN dNims et al.ll2015h . The predicted thermal X-ray luminosity 
at IkeV band is smaller than our non-thermal X-ray prediction and 
dinimishes with increasing outflow shock radius given our assump¬ 
tion about the gas density profile in the galaxy and halo. Since the 
travel time of the outflow shocks is comparable to AGN’s lifetime, 
most of the detected halos still host an active quasar, targets can 
be selected for observations as an AGN. On the other hand, sub¬ 
traction of the much brighter emission from the AGN is required 
to measure the extended diffuse emission from the outflow shocks. 
Radio interferometry can resolve the luminous central source and 
subtract emission from it to obtain the extended emission on halo 
scale. For optical and infrared observations, the extended emission 
can be subtracted usi ng techniques similar to the removal of quasa r 
light in HST images dMcLeod & Riekefl995l : Ba hcall et all 19971) . 

A source of contamination to the extended non-thermal emis¬ 
sion i s the scattered quasar light by the surrounding electrons in th e 
halo dWise & SarazirJ fl 99(ll : I Yound|200(J : iHolder & Loebll2004ll . 
We find that the optical depth for Thomson scattering through 
the halo is ~ ICR 5 , so ~ ICR 5 of the observed flux from the 
AGN is expected to diffuse throughout the halo. For a 10 12 Mq 
mass halo, the bolometric luminosity of the scattered radiation is 
~ 10 40 ergs -1 , which is comparable to the non-thermal emission 
at infrared and optical frequencies from outflow shocks in halos 
within z < 1 and negligible for halos at z ~ 5. One possible way 
to distinguish the scattered radiation from the non-thermal emis¬ 
sion is by polarimetric measurement. Additionally, the scattered 
light is diffused throughout the halo at any given time while the 
emission from outflow shocks shows a discontinuity at the shock 
front. As the outflow propagates farther into the halo, the scattered 
quasar light no longer exists as the quasar fades away. There is 
no contamination from scattered AGN photons in radio band from 
radio-quiet quasars, which takes up ~ 90% of the population, so 
the non-thermal emission can be more easily identified in radio 
wavelength dNims et al.ll2015l) . Therefore, radio observation is ex¬ 
pected to be most effective in detecting the halo scale non-thermal 
emission from outflows in a wide range of redshifts. It should be 
noted that the predicted radio emission from outflow shocks exists 


without the presence of relativ istic jets, which account for the radio 
emission from radio galaxies dHeckman & Bestll2014h . 

There are a few uncertainties in our model. First, spherical 
symmetry of both gas distribution and outflow shell is likely to be 
unrealistic. In fact, the outflow may be collimated from the start or 
can propagate along the path of least resistance, forming a bipo¬ 
lar or bicone structure. Observations of kpc-scale moleculai^out- 
flows suggest a wide-angle biconical geometry dRunke & Veilleu j 
l201ll : iFeruglio et~ai1l2015t) . Biconical outflows with small open¬ 
ing angle could have less impact on the ambient medium. Second, 
the detailed gas distribution is uncertain and can be complicated 
by galaxy-to-galaxy variations, which can greatly dependend on 
galaxy types as well as the specific gas phase. Finally, we find that 
the terminal velocity of the outflow arriving at the edge of the halo 
is ~ 10 3 kms -1 , which is still large enough for farther propaga¬ 
tion of the outflow into the IGM. The propagation dynamics of the 
outflow into the IGM is beyond the scope of this paper. Along some 
directions g as accretion o n to th e galaxy could impede the develop¬ 
ing outflow dSuresh el alJ2015h . 
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